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ABSTRACT 

The  equations  governing  the  motion  of  viscoelastic  materials  with  fading 
memory  incorporate  a  nonlinear  elastic-type  response  with  a  natural 
dissipative  mechanism.  Our  purpose  is  to  discuss  the  subtle  effects  of  this 
mechanism  in  viscoelastic  materials  of  Boltzmann  type.  Recent  results  on  the 
global  existence  and  decay  of  classical  solutions  for  smooth  and  small  data 
(in  one  space  dimension)  are  reviewed  for  smooth  and  singular  memory  kernels; 
for  smooth  kernels  a  number  of  such  results  can  be  generalized  to  several 
space  dimensions.  A  recent  result  on  the  development  of  singularities  in 
finite  time  for  large  data  is  discussed;  several  open  problems  are 
formulated.  A  program  for  studying  weak  solutions  for  such  systems,  including 
the  development  of  numerical  algorithms,  is  outlined.  -  • 
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NONLINEAR  VISCOELASTIC  MATERIALS  WITH  FADING  MEMORY 


John  A.  Nohel 

1.  Introduction.  The  equations  governing  the  motion  of  nonlinear  elastic  bodies  are 
quasilinear  hyperbolic  systems  for  which  smooth  solutions  generally  lose  regularity  in 
finite  time  due  to  the  formation  of  shock  fronts.  Some  materials  incorporate  a  nonlinear 
elastic-type  response  with  a  natural  dissipative  mechanism,  and  it  is  important  to 
understand  the  effects  of  the  dissipation  on  the  behaviour  of  the  solutions  of  the 
equations  of  motion. 

The  purpose  of  this  lecture  is  to  discuss  the  effects  of  the  subtle  dissipative 
mechanism  due  to  memory  effects  in  viscoelastic  materials  of  Boltzmann  type.  This 
dissipation  is  more  delicate  than  that  exhibited  by  viscoelastic  materials  of  the  rate  type 
for  which  globally  defined  smooth  solutions  exist,  even  for  large  smooth  data. 

The  paper  is  organised  as  follows.  In  Section  2  we  formulate  mathematical  models  for 
the  notion  of  nonlinear  viscoelastic  materials  and  we  motivate  the  mathematical  theory.  In 
Section  3  we  survey  recent  results  on  the  global  existence  of  smooth  solutions  for  smooth 
and  small  data.  In  Section  4  we  present  a  recent  result  on  the  breakdown  of  smooth 
solutions  for  large,  smooth  data  and  discuss  briefly  related  open  queetions  including  those 
regarding  weak  solutions  and  numerical  methods  (Remarks  4.8).  We  restrict  our  attention 
throughout  to  one-dimensional  problems  and  provide  some  referencee  for  multidimensional 
problems.  Moreover,  we  consider  only  a  purely  amchanlcal  theory,  i.e.  we  neglect  thermal 
effects. 


2.  Mathematical  Models  and  Dynamic  Problems.  Consider  the  longitudinal  motion  of  a 
homogeneous  one-dimensional  body  (e.g.  a  bar  of  uniform  cross-section)  occupying  an 
interval  B  in  a  reference  configuration,  which  we  assume  to  be  an  equilibrium  state,  and 
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having  unit  reference  density*  B  aay  be  bounded  or  unbounded.  Let  u(x,t)  denote  the 
displacement  at  time  t  of  a  particle  with  reference  position  x  (l.e.  x  +  u(x,t)  Is  the 
position  at  time  t  of  the  particle  at  x).  The  strain  which  measures  local  stretching  Is 


defined  by  c  =  ux(x,t).  Let  o  denote  the  stress  at  time  t  of  the  particle  with 
reference  position  x  (o  Manures  the  contact  force  per  unit  area).  The  balance  of  linear 
momentum  yields  the  equation  of  motion 

utt  “  °x  *  f  '  x  «  B'  >  0  * 

where  subscripts  denote  partial  derivatives  and  where  f  is  an  external  body  force.  In 
order  to  characterise  the  material,  (2.1)  is  supplemented  by  a  constitutive  assumption 
which  relates  the  stress  to  the  motion.  In  addition,  initial  data,  as  well  as  suitable 
boundary  data  if  B  is  not  R,  are  adjoined  to  (2.1).  We  remark  that  in  a  physical 
problem  the  cross-section  does  not  generally  remain  uniform  as  the  bar  is  stretched.  More 
realistic  problems  can  be  treated  by  similar  techniques. 

If  the  body  is  homogeneous  and  purely  elastic,  the  stress  depends  on  the  strain 
through  the  constitutive  relation  o(x,t)  »4(e(x,t)),  where  4  is  a  given  smooth 
function  satisfying  the  assumptions  (i)  4(0)  -  0,  (il)  ♦’(0)  >  0»  (i)  reflects  the  fact 
that  the  reference  position  is  taken  as  an  equilibrium  state,  and  (ii)  that  the  stress 
increases  with  the  strain,  at  least  near  equilibrium.  The  equation  of  motion  (2.1)  becomes 
the  familiar,  one-dimensional,  quasllinear  wave  equation 

“tt  “  ♦(ux,x  *  f  (x  «  »,  t  >  0)  (  (2.2) 

if  B  is  bounded  it  is  assuMd  that  the  assigned  boundary  data  and  initial  data  are 
compatible.  For  (2.2)  there  is  no  natural  dissipative  Mchanlsm.  Indeed,  Lax  (33],  also 
MacCamy  and  Nisei  (37]  and  Kleinerman  and  Hajda  [31]  have  shown  that  if  4  is  not  linear, 
the  Cauchy  problem  for  (2.2)  (f  *  0)  does  not  generally  possess  globally  defined  smooth 
solutions,  no  matter  how  smooth  and  saall  one  takes  the  initial  data  u(x,0)  and  u^Cx.O). 

In  a  material  with  memory  (such  as  certain  polyMrs,  suspensions,  or  emulsions)  the 
stress  at  a  material  point  x  and  at  time  t  depends  on  the  entire  history  of  the  strain 
at  x.  In  1874  Boltzmann  [5]  gave  the  following  linear  constitutive  law  for  saall 


m 
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deformations  in  such  materials 

c(x,t)  -  Be(x,t)  +  /”  m(s) [e (x,t)  -  e (x,t-s)]ds,  xeB,  -  «  <  t  <  “  .  (2.3) 

In  (2.3)  8  >  0  is  a  given  constant  and  m  i  (0,“)  ♦  R  is  a  given  positive,  smooth, 
nonincreasing  function.  We  limit  our  discussion  to  the  situation  in  which  m  c  L*(0/»), 
and  we  distinguish  two  cases: 

(i)  0  <  m( 0)  <  «,  (ii)  m(0+)  -  .  (2.4) 

The  function  m  is  called  a  memory  function.  The  fact  that  m  >  0  and  non-increasing 
on  (0,»)  means  that  the  stress  'relaxes"  as  t  increases  and  the  memory  term  in  (2.3) 
fades:  deformations  which  occurred  in  the  distant  past  have  leas  Influence  on  the  present 
value  of  the  stress  than  those  which  occurred  in  the  recent  past.  In  the  applied 
literature  m  is  often  assumed  to  be  a  finite  linear  combination  of  decaying  exponentials 
with  positive  coefficients  (these  expressions  result  from  least  squares  approximations  to 
experimental  data).  Such  restrictions  are  neither  desirable  nor  necessary.  Moreover, 
kinetic  theories  for  chain  molecules  [15,46,53]  and  certain  experiments  [32,28]  suggest 
that  there  are  materials  for  which  m  is  singular  as  in  (2.4)(li),  m(t)  ~  ta-1  as 
t  +  0+,  0  <  o  <  1,  m  Is  positive,  nonincreasing  on  0  <  t  <  •*,  and  m  decays  rapidly  at 
infinity.  Stronger  power  singularities  at  zero  (a  <  0)  are  also  possible,  but  the 
resulting  mathematical  theory  for  nonlinear  materials  consistent  with  our  objectives  is 
incomplete  at  this  time. 

The  assumption  m  e  L^O,*)  implies  that  (2.3)  is  equivalent  to 

o(x,t)  -  c2e(x,t)  -  /p  m(s)e  (x,t-s)da,  xeB,  -  »  <  t  <  »  ,  (2.5) 

where  c2  :»  B  +  /”  m(s)ds  >0  is  a  constant  which  measures  the  instantaneous  response  of 
stress  to  strain:  8  >  0  is  the  equilibrium  stress  modulus.  If  0  >  0  the  material  acts 
like  a  solid,  while  if  8  ■  0  it  acts  like  a  fluid. 

A  natural  generalization  of  (2.5)  to  nonlinear  materials  is  the  constitutive  relation 
o(x,t)  -  4(e(x,t))  -  m(s)*(e(x,t-s))ds,  xeB,  —  <  t  <•  ,  (2.6) 

in  which  ♦ ,  *  :  R  +  R  are  assigned,  smooth  material  functions  which  satisfy 

♦  (0)  -  *(0>  -  0,  ♦'(0)  >  0,  i|>’<0)  >  0  .  (2.7) 
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The  memory  function  ra  ia  positive,  nonincreasing  and  integrable  on  (0,“)  as  above.  In 
the  static  case  e(x,t)  «  c(x),  o(x,t)  »  a(x),  (2.6)  reduces  to 

a(x)  -  ♦  (e  (x) )  -  </*  m(  s Ids  )4>  («T(x) ) ,  x  c  B  . 

A  natural  assumption,  appropriate  for  viscoelastic  solids  and  c.uclal  in  the  analysis  of 
global  existence  results  (section  3),  is  to  require  that  + ,  also  satisfy 

♦  ’<0)  -  (/q  m(s  )ds)ii '  (0)  >  0  *  (2.8) 

(2.8)  states  that  the  equilibrium  stress  modulus  is  positive.  The  constitutive  assumption 
(2.6)  is  a  particular  case  of  a  "simple  material"  (8)  which  retains  many  important 
qualitative  properties  of  more  general  material  models*  moreover,  the  analysis  of  the 
resulting  equation  of  motion  is  relatively  simple  and  complete. 

The  balance  of  linear  momentum  and  (2.6)  yield  the  equation  of  motion 

utt  *  ♦‘Vx  "  /l o  »<t-T)*(ux(x,T)>x4r  +  f,  xeB,  —  <t<»  ,  (2.9) 

where  f  is  a  body  force  and  where  the  change  of  variable  t  :»  t-s  was  made  in  (2.6). 

The  history  of  the  motion  is  assumed  to  be  known  for  t  <  0  (the  history  may,  but  need  not 
satisfy  (2.9)  for  t  <  0).  An  appropriate  dynamic  problem  is  to  find  a  smooth  function 
u  s  B  x  (-“,“)  ♦  R,  satisfying  (2.9)  for  t  >  0,  and  such  that 

u(x,t)  ”  u(x,t)  ,  x  e  B,  t  <  0  ,  (2.10) 

where  the  history  u  s  B  *  (-“,0]  R  is  a  given  smooth  function*  (2.9),  (2.10)  will  be 
referred  to  as  a  history  value  problem.  If  B  is  bounded  or  semibounded  compatible 
boundary  conditions  are  adjoined  to  (2.9),  (2.10).  Compatibility  of  the  boundary 
conditions  with  the  smooth  data  f  and  u  is  imposed  in  order  to  preclude  the  propagation 
of  singularities  from  the  boundary  into  the  interior. 

If  ra  3  0,  (2.9)  reduces  to  the  quasilinear  wave  equation  (2.2).  At  the  other 
ixtreme,  if  one  formally  sets  m  «  -6  * ,  where  6  is  the  Dirac  mass  at  the  origin,  then 

(2.9)  reduces  to  the  parabolic  equation 

utt  "  *<ux>xt  +  ♦'Vx  +  f  > 

the  terra  ♦(ux)xt  represents  viscosity  of  Newtonian  type  if  ip  is  smooth  and  *p  ’  (•  )  >  0. 


This  equation  possesses  globally  defined  smooth  solutions  even  if  the  data  are  large 


(1.34J. 

Our  objective  is  to  discuss  the  strength  of  the  dissipative  mechanism  in'*  iced  by  the 
memory  in  (2.9)  under  physically  reasonable  assumptions  by  studying  the  existence  and  the 
decay  or  growth  of  classical  solutions  of  the  history  value  problem  (2.9),  (2.10).  To 
motivate  the  mathematical  results,  we  follow  Coleman  and  Gurtin  [6]  in  their  penetrating 
study  on  the  growth  and  decay  of  acceleration  waves  propagating  into  a  one-dimensional 
viscoelastic  material  with  memory  at  rest.  An  acceleration  wave  solution  u  is  similar  to 
a  shock  wave;  the  difference  is  that  second  rather  than  first  derivatives  of  u  experience 
a  jump  acrosss  the  wave  front.  To  apply  the  results  of  [6]  to  (2.9),  (2.10),  we  assume 
that  4,  4>  are  smooth,  satisfy  (2.7),  f  S  0,  B  -  R,  and  m  is  a  smooth,  regular  kernel 
satisfying  (2.4)(i).  “Hie  wave  front  is  a  smooth  curve  t  «  Y(x),  7(0)  «  0,  and  u  =  0 
for  t  <  y(x).  In  [6]  the  problem  of  existence  of  acceleration  waves  is  not  discussed. 
Assuming  that  they  do,  an  easy  but  tedious  calculation  shows  that  for  (2.9)  t  -  y (x)  is  a 
straight  line,  of  slope  (♦ * (0) )“  ^ ,  meaning  that  such  waves  propagate  with  constant 
speed  although  (2.9)  is  nonlinear.  I<et  the  amplitude  of  the  wave  be  q(t)  [uttl , 
where  (ufct]  is  the  jump  in  across  the  line  t  “  y(x).  It  follows  from  the 

computations  in  (6)  that  q  evolves  in  accordance  with  the  Ricatti-Bernoulll  equation 

ft  *  -  Aq2  -  Bq  ,  q(0)  -  qQ  ,  (2.11) 
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where  —  ”  r —  ♦  c  -r—,  c  ”  4’{0),  represents  differentiation  along  the  wave  front  and 
dt  it  3x 
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given  by  [ 

solution  of  (2.11)  tends  to  zero  as  t  +  +».  By  contrast,  if  q0  >  +  B/A,  then  q(t)  ♦  ** 


^0 

»VB 


>  0.  the  corresponding  jumps  in  xi%t  and  uxx  are 


1  W 

given  by  [uxt)  “  - M 0) ] ^ q( t)  and  [u^]  “  [♦'(0)]"qt). 


This  rssu' t  suggests  the  following  conjectures  regarding  smooth  solutions  of  the 


history  value  problem  (2.9),  (2.10): 

(i)  The  problem  (2.9),  (2.10)  should  »-  -we  globally  defined  classical  solutions  if  the 
history  v  and  the  forcing  term  f  are  sufficiently  smooth  and  small  in  appropriate 
norms.  Moreover,  such  solutions  should  decay. 

(ii)  The  smooth  solutions  of  (2.9),  (2.10)  should  develop  singularities  in  second 
derivatives  in  finite  time  if  the  smooth  data  are  chosen  sufficiently  large. 

As  will  be  summarized  in  Section  3,  conjecture  (i)  has  been  established  rigorously  by 
a  number  of  authors  in  a  number  of  physically  important  cases  of  (2.9),  (2.10)  for  regular 
kernels  (m(0)  <  »),  as  well  as  for  singular  kernels  (m(0+)  »  ■*<•).  Conjecture  (ii)  has 
only  been  established  for  regular  kernels  (see  Section  4).  Moreover,  based  on  the 
discussion  in  Section  4,  Remark  4.6,  singular  kernels  m  strengthen  the  dissipative 
mechanism  of  the  memory  in  (2.9)  which  suggests  the  possibility  that  for  appropriate 
classes  of  singular  kernels,  global  smooth  solutions  will  exist  even  if  the  data  are 
arbitrarily  large i  this  interesting  question  is  open. 

Moat  of  the  results  described  in  Sections  3  and  4  for  smooth  kernels  satisfying  (2.4a) 
apply  to  more  general  one-dimensional  viscoelastic  models  with  fading  memory,  e.g.  a  model 
for  a  solid,  K-BKZ  material  [29,2] 

u..  -  4(u  )  +  ra(t-T)h(u  (x,t),u  (x,t  ) )  dr  ♦  f,  (2.12) 

tr  A  ^  x  xx 

xeB,  -*»<t<«* 

Here  4,  m,  and  f  are  as  in  (2.9),  while  h  :  R  x  R  ♦  R  is  a  smooth  material 
function,  h(p,p)  «  0  and  the  partial  derivatives  of  h  satisfy  appropriate  sign 
conditions,  at  least  at  (0,0).  If  4=0,  (2.12)  models  a  K-BKZ  fluid.  Under  suitable 

assumptions,  the  energy  method  for  proving  existence  results  in  Section  3  and  the  method  of 
characteristics  used  to  prove  blow-up  results  of  Section  4  yield  similar  results  for  this 
case  as  well.  The  energy  method  can  also  be  applied  to  prove  existence  for  certain 
multidimensional  viscoelastic  problems  with  fading  memory  (e.g.  [13,  Sec.  4],  [30]). 
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However,  to  our  knowledge,  the  existence  results  described  in  Section  3  for  singular 
kernels  satisfying  (2.4(ii))  depend  crucially  on  the  special  fora  of  equation  (2.9). 


I 

i 
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3.  Existence  of  Classical  Solutions.  For  discussion  of  the  mathematical  results  it 
is  convenient  to  renormalize  the  memory  function  m.  Define  the  relaxation  function  a  by 

a(t)  i-  /“  m ( s ) ds ,  0  <  t  <  “  i  (3.1) 

observe  that  if  m  is  smooth,  positive,  decreasing  and  lntegrable  on  [0,“)  then 
a' (t)  -  -m(t)  and 

a  is  smooth,  positive,  decreasing  and  convex  on  (0,»)  .  (3.2) 

Analogous  to  (2.4)  we  distinguish  two  classes  of  kernels  as 

(i)  0  <  -a' (0+)  <  -  ,  (ii)  -a * ( 0+ )  -  t-  .  (3.3) 

other  normalizations  of  the  memory  m  are  possibles  for  example,  the  relaxation  function 
G(t)  4'(0)  -  s(0)*'(0)  +  a(t)**(0)  ,  0  <  t  <  -  ,  (3.4) 

where  t ,  <|>  are  the  material  functions  in  (2.6),  is  consistent  with  the  applied 
literature.  Observe  that  G(“)  -  ♦'(<))  -  a(0)(i'(0)  and  G(0)  «4'(0). 

Returning  to  the  history  value  problem  (2.9),  (2.10),  let  the  history  u  be 
identically  zero  for  t  <  0.  One  then  seeks  a  solution  of  the  initial  value  problem 

utt  -  ♦<uJC)x  +  /g  aMt-r)*(tt  (x.T))^*  +  f,  x  €  B,  t  >  0  ,  (3.5) 

u(x,0)  -  Uq(x),  Uj.tx.O)  »  u,(x),  X  c  B  ,  (3.6) 

together  with  suitable  and  conqsatible  boundary  conditions  if  B  is  not  R.  If  the  history 
u  is  not  zero  for  t  <  0,  the  part  of  the  Integral  in  (2.9)  on  (-",0)  is  Incorporated 


Global  Existence  of  Classical  Solutions.  He  next  discuss  global  existence  and 
asymptotic  behaviour  for  the  Cauchy  problem  (3.5),  (3.6)  with  B  «  R,  for  smooth,  small 
data,  and  for  regular  kernels  a  satisfying  (3.2),  (3. J)(i).  To  simplify  the  exposition, 
we  make  the  hypothesis 


•‘%UV 


a  e  crto,-),  (-1)  a'K'(t)  >0(0<t<-ik-  0, 1,2,3)  , 

a*  f  0,  and  ta(t)dt  <  ® 


The  results  hold  under  assumptions  on  a  considerably  weaker  than  (3.7).  The  interested 
reader  is  referred  to  [131,  (221,  [24],  and  the  survey  paper  [23]  for  the 
generalizations.  The  essential  point  is  that  kernel  a  satisfies  a,  a',  a"  c  L^O,-), 
the  moment  condition  in  (3.7),  and  is  "strongly  positive"  on  [0,“).  The  result  for  B 
bounded  [13]  is  somewhat  simpler  than  for  the  Cauchy  problem  (3.5),  (3.6);  in  particular, 
the  moment  condition  in  (3.7)  is  only  needed  for  the  Cauchy  problem  (see  remarks  following 
Theorem  3.1  and  the  outline  of  its  proof). 

Concerning  $ ,  i(i  assume 

♦  ,  <l>  £  C3(«),  4(0)  -  *(0)  -  0,  ♦•(0)  >  0,  *  •  ( 0)  >  0,  ♦•(0)  -  a(0)$  '  (0)  >  0  t  (3.8) 

the  latter  is  the  analogue  of  (2.8)  in  the  present  normalization.  Assume  that 


(i)  f,  fx,  ffc  £  C([0,«)»  L2(B))  D  l“([0,-)i  L2(*>)  and 

(ii)  f  £  L1([0,-)»  L2<«)),  f  ,  f  ,  f  c  L2([0,-)»  L2(*)) 

X  t  xt 


(3.9) 


and  let  uQ,  Uj  satisfy 

U0  €  Ltoc(K)'  *nd  U0  U1  C  h2(R>  *  (3-10) 

To  measure  the  size  of  the  data  define  the  quantities 

U0<W  /_„{Ug2  ♦  UJ2  +  u5'2  +  u2  ♦  u;2  +  u"2}  (x)dx,  and  (3.11) 

P(f)  sup  f2  t  f2  +  f2}(x,t)dx  +  (/"(/^  f2(x,t)dx)1/2dt)2  (3.12) 

t>0  x  t  Q  — 

+  /q  U'l  *  ft  +  ^t}(X't)dXdt  ' 


The  following  result  is  a  special  case  of  Theorem  1.1  of  [24]. 


Theorem  3.1.  Let  assumptions  (3.7)  -  (3.10)  be  satisfies.  There  exists  a  constant 
U  >  0  such  tnat  for  each  u0,  uj,  f  satisfying 

U(u0,u,)  +  F(f)  <  u2  ,  (3.13) 


-8- 


the  Cauchy  problem  (3.5),  (3.6)  has  a  unique  solution  u  €  C^(E*  1 0  /»  ) ) ,  and 


Moreover , 


u  ,  u  ,  a  - - ,u  £  C([0,«);  L2(R))  n  l"([0,»)j  L2(K)) 

X  t  XX  ttt 


uxx'  uxt . uttt  C  L2([0,-),  L2(«)  , 


(3.14) 


(3.15) 


u  ,  u  ,  u,_^  *  0  in  L  (*)  as  t 
XX  xt  tt  — 


(3.16) 


u  ,  u.  ,  u  ,  u  .  ,  u.  .  ♦  0  uniformly  on  R  as  t  ♦  ■ 

X  t  XX  xt  tt  - * -  — 


(3. 17) 


A  similar  result  holds  for  the  history  value  problem  (2.9),  (2.10)  with  B  =  R.  The 
special  case  a(t)  =  ae-^,  a  >  0,  X  >  0,  studied  by  Greenberg  [18]  for  B  bounded,  is 
carried  out  in  [23]  in  the  more  complicated  case  when  B  =  R. 

Remark  3.2.  Theorem  3.1  is  a  generalization  of  Theorems  1.1  and  4.1  of  [13]  establishing 
small-data  global  existence  results  for  analogous  initial  boundary  value  problems 
corresponding  to  motions  of  bounded  viscoelastic  bodies;  Neumann,  Dirichlet  and  mixed 
boundary  conditions  are  treated.  The  principal  difficulty  in  proving  Theorem  3.1  is  that 
various  Poincare  inequalities,  not  applicable  to  (3.5),  (3.6)  when  B  =  R,  are  used  in 
an  essential  way  in  [13]  to  establish  an  a  priori  estimate  similar  to  (3.26)  from  (3.32) 
(see  outline  of  proof  following  Proposition  3.4);  the  estimate  (3.26)  is  essential  for 
completing  the  proof.  The  reader  is  referred  to  Hrusa  [22]  for  a  discussion  of  general 
history  value  problems  on  a  bounded  interval.  Although  technically  extremely  complicated, 
the  generalization  of  the  results  in  [22]  to  the  Cauchy  problem  is  relatively  straight¬ 
forward. 


Remark  3.3.  If  i|i  =  ♦  equations  of  the  form  (3.5)  have  been  studied  by  MacCamy  [35J, 
Dafermos  and  the  author  [12],  and  Staffans  [49]  for  bounded  and  unbounded  bodies.  If 


=  p,  (3.5)  admits  certain  estimates  which  do  not  carry  over  to  the  general  case  i|>  ?  * 
(see  [23]);  there  does  not  appear  to  be  any  physical  motivation  for  the  restriction  ♦  =  ♦ 
for  solids. 


Outline  of  the  proof  of  Theorem  3.1.  An  essential  ingredient  of  any  global  result  is  an 


appropriate  local  existence  theorem.  For  regular  kernels  a  satisfying  (  3.  2 ) ,  (  3.  3)  ( i ) , 


the  idea  is  to  iterate  the  sequence  of  linear  problems  which  treat  the  memory  as  a  lower- 


order  perturbation: 


u.  .  -  *  ’  (w  )u„v  +  /5  a’  (t-T  )i|i(w  (x,T  ) )  dr  +  f ,  xe  R,  0  <  t<  T  ,  (3.18) 

tt  A  U  XX 


where  T  >  0,  u  satisfies  the  initial  conditions  (3.6),  and  where  w  is  an  element  of  a 


suitably  chosen  function  space  X.  By  using  fairly  standard  enery  estimates  deduced  from 


(3.18),  requiring  only  very  simple  estimates  of  the  convolution  term  which  do  not  use  any 


sign  information  on  the  memory,  it  is  shown  that  the  mapping  S  which  carries  w  into  a 


solution  of  (3.18)  has  a  unique  fixed  point  for  T  >  0  sufficiently  small.  The  proof  is 


almost  identical  with  that  of  Theorem  2.1  of  [13].  The  only  significant  difference  is  that 


the  proof  in  [13]  is  for  x  €  [0,1]  with  Neumann  boundary  conditions  satisfied  at  x  «  0 


and  x  =  1 )  thus  the  Poincare  inequality  enables  one  to  deduce  estimates  for  lower  order 


derivatives  of  u  in  L°°([0,T)|  L2(0,1))  from  higher  order  derivative  estimates.  As  far 


as  local  existence  is  concerned  when  B  =  R,  this  causes  no  serious  difficulties.  One 


simply  expresses  the  lower  order  derivatives  of  the  solution  in  terms  of  initial  conditions 


and  time  integrals  of  the  higher  order  derivatives,  yielding  time  dependent  bounds  which. 


however,  cannot  be  used  for  obtaining  global  estimates.  The  result  is: 


Proposition  3.4.  Let  a,  a',  a"  c  L^^IO,**)  and  assume  that 


£  CJ(E), 


>  •  ( 0 )  >  0,  and  that  there  exists  a  number  such  that 


*'(5)  >  ♦  for  every  ?  e  R  • 


(3.19) 


Concerning  the  data,  let  Ug,  u1  satisfy  (3.10),  f  satisfy  (3.9)(i)  and 


1  2 

assume  that  fxt  €  Lf0c^ 0 >'’’)!  L  (R)).  Then  the  Cauchy  problem  (3.5),  (3.6)  has  a  unique 


solution  u  defined  on  a  maximal  time  interval  [ 0 , Tn )  satisfying 


V  V  Uxx'  Uxt'  utt'  uxxx'  uxxt'  \tt'  Uttt  £  C('°'V’  L  (R))  •  (3-20) 


Moreover ,  i f 


t£(0,T0) 


O-J  +  Ut  +’”  +  uxtt  +  uttt)(x't)dx  <  *  ' 


(3.21) 


•gw 

*!w 


then  Tn  «  +“*.  By  Sobolev  embedding  u  e  C  (R  x  [0,T  )). 


In  outline,  the  proof  of  the  global  result  then  proceeds  as  follows.  Define  the 
equilibrium  stress  x  by 


X(€)  s-  ♦  (£)  -  a(O)iHC),  t  I  i 


(3.22) 


observe  that  x  £  C  (R)  and  that  x*(0)  >  0  (by  3.8).  Choose  a  sufficiently  small 
number  5  >  0  and  modify  ♦  ,  ii ,  and  X  outside  [-6,5]  such  that  i|)“,  X"  vanish 

outside  [-25,25],  and  choose  positive  constants  £,  ^  such  that 


*’(5)  >  1,  *'(€)  >  £,  x'(5)  >  1  e  R 


(3.23) 


It  is  shown  a  posteriori  that  |ux(x,t)|  <  5  for  all  x  e  R,  t  >  0.  By  Proposition  3.4 
the  Cauchy  problem  (3.5),  (3.6),  B  ■»  R  has  a  unique  solution  u  on  a  maximal  interval 
[ 0,Tq ) .  The  objective  is  to  show  that  if  (3.13)  holds  with  y  >  0  sufficiently  small, 
then  (3.21)  is  bounded  independent  of  Tg»  a  standard  continuation  procedure  implies 


Tq  =  +<■>.  Define 


E(t)  max  /"  {u2  +  u2  +•••+  u2  )  (x,s)dx 

s€[0,t]  *  * 

+  ^0  {vi  *  +*,,+  uLt}(x'a)dxd»  ' 


(3.24) 


where  • • •  represent  the  sum  of  the  second  and  third  derivatives  not  explicitly  written 
down.  It  is  shown  that  if  (3.13)  holds  for  y  >  0  sufficiently  small,  then  E(t)  is 


bounded.  For  this  purpose  define 


,2  *  2  *  2  iV2 


v(t)  sup  {ux  +  uxx  +  u^)  z(x,s),  Vt  c  [0,T0)  . 

xe  R  J 

s<  ( 0 ,  t  ] 

To  prove  the  result  one  establishes  the  following  Key  estimate 

E(t)  <  r(U0(uQ,u1)  +  F(  f ) )  +  T(v(t)  ♦  v3(t)  )E(  t ) ,  0  <  t  <  TQ  , 


(3.25) 


(3.26) 


where  here  and  below  r  is  a  generic  constant,  possibly  large,  independent  of  uQ,  u^> 
f,  and  Tg.  We  shall  comment  below  only  briefly  how  this  is  accomplished. 

Once  (3.26)  is  established,  the  conclusions  of  Theorem  3.1  are  obtained  as  follows. 


Choose  E,  y  >  0  such  that 


e  <  52,  r{(2E)1'^2  +  (2E)3/2}  <  7,  ry2  <  jE  . 


Select  the  data  Ug ,  Uj,  f  such  that  (3.13)  holds  for  u  chosen  in  accordance  with 


rot 

S'? 

PS 

l  LT-  ‘2 


(3.27).  The  Sobolev  embedding  theorem  implies  that 


v(t)  <  (2E(t))/2  Vt  e  t°*T0>  • 


(3.28) 


Therefore,  it  follows  from  (3.26),  (3.27),  (3.28)  that  for  any  t  r  [0,Tg)  with 

E(t)  <  E,  we  actually  have  E(t)  <  j  E.  By  continuity  E(t)  <  y  E,  »t  £  [ 0 , Tg ) , 

1  —  o 

provided  E(0)  <  —  E)  the  latter  is  insured  by  choosing  smaller  if  necessary  so  that 

(3.13)  will  imply  E(0)  <  -y  E.  Then  E(t)  <  y  E ,  Vt  £  [0,Tg),  and  (3.24),  Proposition 
3.1,  and  a  standard  continuation  method  yield  Tg  «  +“.  One  also  has  that  (3.14),  (3.15) 
hold,  and  conclusions  (3.16),  (3.17)  follow  by  standard  embedding  inequalities.  Moreover, 
(  3.25),  (3.27),  (  3.28)  yield 

|ux(x,t)|  <  v(t)  <  (2E(t))1/2  <  (I)1/2  <  6,  Vx  £  R,  t  e  (0,-) 
and  the  proof  is  complete. 

Establishing  the  estimate  (3.26)  is  lengthy,  delicate,  and  relies  on  the  correct  sign 
of  the  memory  (under  assumption  (3.7)  or  certain  generalizations].  The  energy  method, 
combined  with  relevant  properties  of  Volterra  operators  and  their  resolvents,  is 
employed.  The  estimates  of  derivatives  of  u  appearing  in  (3.24)  are  deduced  from  energy 
identities  obtained  directly  from  (3.5),  (3.6),  and  from  the  equation  equivalent  to  (3.5): 


|r>  :• 


Utt  -  X<Vx  +  /;  »(t-T)*(Ux)jit(x,T)dT 


(3.29) 


+  a(t)*(u0(x)x)x  ♦  f  (x  e  *,  0  <  t  <  T)  , 


where  T  <  Tg)  (3.29)  is  obtained  from  (3.5),  (3.6)  by  an  integration  by  parts  and  use  of 
(3.22).  Useful  Identities  for  derivatives  of  u  can  only  be  obtained  by  multiplying  the 
equations  by  quantities  which  make  it  possible  to  estimate  the  memory  terms,  h  crucial 
role  is  played  by  the  "quadratic  integral  form" 

Q(w,t,b)  :»  /g  FI*  w(x,s)  J*  b(s-x  )w(x,t  )drdxds,  t  >  0  , 

1  2 

defined  for  b  e  an<*  ^or  avery  w  e  c(  (0»t)  )  L  (*)).  In  the  first  energy 


and  b  »  a.  It  is  an  important  fact  that 


over  R  x  (0,T]  ,  Q  arises  with  w  =  ♦  (Ujj)^ 


kernels  a  satisfying  (3.7)  (indeed  much  weaker  assumptions)  are  positive  definite  on 
(0,®).  To  obtain  the  second  energy  identity,  one  needs  to  take  the  forward  time-difference 
of  (3.29)  and  integrate  the  resulting  equation  over  R*  [0,T]  .  To  estimate  the  relevant 


derivatives  of  u  from  a  combination  of  the  first  two  identities  one  needs  the  following 
technical  estimate:  It  is  shown  in  [24i  Lemma  2.5]  that  if  a  satisfies  (3.7),  there 
exists  a  constant  x  >  0  such  that 


f„  /"  w2(x,t)dxdt  <  x  /  w  (x,0)dx  +  xQ(w ,t,a) 
0  *  — t  ■<*  t  t 


(3.30) 


+  x  lim  inf  —  C(A  w  ,t,a)  ,  V  t  e  [0,T]  , 

,  ,  _  ,  Z  rl  t 

h+0  h 


1  2 

where  i  C  ( [0,T] »  L  (E))  V  T  >  0,  and  where  the  forward  difference  operator  Ahw  is 
defined  by  Ahw(x,t)  :«  w(x,t+h)  -  w(x,t).  In  the  application  of  (3.30),  w  ■>  i|j  (ux  *xt  and 
the  forward  difference  operator  A^  is  applied  to  equations  (3.29).  The  proof  of  (3.30) 
also  makes  use  of  a  result  of  Staffans  ({49,  Lemma  4.2]).  Using  the  two  energy 
identities,  and  (3.30),  it  is  relatively  straightforward  to  estimate  all  of  the  terms  and 


arrive  at: 


f“  (u2  ♦  u2  +  u2  +  u2  }(x,t)dx  +  /”  u2  (x,s)dxds 

‘  -»  xx  xt  xxt  xtt  ‘  0  '  -»  xxt 

«  r(u0+-n  +  r(v(t)+v3(t) >s(t)  +  r(/u^  +  /?)  /iiftT,  w  t  t  [o,t]  . 


(3.31) 


Estimates  of  u2t(x,t)dx,  u2fct(x,t)dx,  /  —  u2tfc(x,T  )dxdt,  V  t  e  [0,T]  in  terms 

of  the  right  side  of  (3.31)  are  obtained  from  (3.5).  A  bound  for  u2t(_(x,s)dxds 

can  then  be  obtained  by  interpolation.  Using  the  fact  that  a  certain  resolvent  kernel  of 
a'  in  (3.5)  is  in  L’tO,®),  Lemma  3.2  of  [13]  makes  it  possible  to  estimate 
f  u2  (x,t)dx  and  f!r  /"*  u2  (x,s)dxds.  Combining  these  with  (3.31)  yields  the 

estimate 


m. 


I  lu  +  u  +  u  ♦  u  +  u  .  +  u  +  u.  >  (x,t)dx  (3-32) 

J  -•  XX  Xt  tt  XXX  xxt  xtt  ttt 

+•  f  {u2  +  u2  +  u2  +  u2  }  (x,a)ds 

J  0  J  —  XXX  xxt  xtt  ttt'  ' 

<  r(uQ  +  f>  +  r(v(t)  +  v3(t))E(t) 

+  r(/U^  +  /F)  /E(t),  V  t  £  [0,T] 

The  estimate  (3.32)  is  implicit  in  the  argument  of  [13].  It  should  be  observed  that 
for  problems  on  bounded  intervals  (Remark  3.2),  it  is  a  simple  matter  to  apply 
the  Poincare  inequality  to  deduce  the  remaining  estimates  of  derivatives  of  u  appearing 
in  (3.24)  and  arrive  at  the  final  estimate  (3.26)  directly  from  (3.32).  However,  to 
accomplish  this  task  for  (3.5),  (3.6)  when  B  »  R  is  quite  tricky  and  involves  additional 
properties  of  Volterra  operators  and  certain  other  of  their  resolvents.  The  reader  is 
referred  to  Lemmas  2.3  and  2.4,  as  well  as  the  argument  on  pages  405-410  of  [24]  for 
details.  This  part  of  the  proof  makes  essential  use  of  the  assumption  a"  £  L 1 ( 0 ,"  )  which 
is  automatic  when  a  satisfies  (3.7),  txit  cannot  be  satisfied  by  singular  kernels. 

For  singular  kernels  satisfying  (3.2)  and  (3.3)(ii),  it  is  simpler  to  restrict  the 
analysis  to  the  history  value  problem,  (2.9),  (2.10),  with  a  defined  by  (3.1),  in  which 
that  history  u  satisfies  the  equation  (and  the  boundary  conditions  if  B  is  bounded). 
This  ensures  that  the  compatibility  conditions  between  the  history  and  boundary  data,  as 
well  as  compatibility  conditions  between  the  derivatives  of  the  history  and  the  solution 
for  t  >  0  are  satisfied.  If  u  is  a  smooth  solution  of  (2.9)  and  the  kernel  a  is 
singular,  the  Integral  in  (2.9)  is  also  a  smooth  function,  but  the  integrals  /5L  and  /q 
have  singularities  at  t  «  0  which  cancel.  Thus  if  formulated  as  an  initial  value  problem 
the  results  would  involve  a  singular  forcing  term.  For  reasons  explained  below,  global 
existence  results  for  singular  kernels  only  hold  for  B  bounded. 

The  principal  difficulty  when  dealing  with  singular  kernels  is  establishing  a  suitable 
local  existence  result.  In  Proposition  3.4  for  regular  kernels  no  hypothesis  is  made 
concerning  the  sign  of  the  memory  and  the  size  of  the  data.  In  the  proof  the  memory  is 
treated  as  a  perturbation  of  the  elastic  terra  ^  ( va  )  in  (3.5).  However,  the  proof  makes 


a 


crucial  use  of  the  hypothesis  a"  e  which  rules  out  singular  kernels  a 

satisfying  (3.2),  (3.3)(ii). 

Hrusa  and  Penardy  [25,  <n\eorem  4.1]  recently  obtained  an  elegant  extension  of 
Proposition  3.4  for  such  singular  kernels.  They  consider  the  history  value  problem  with 
the  history  satisfying  the  equation  and  the  boundary  conditions  for  t  <  0.  The  singular 
kernel  a  satisfies  the  assumptions 

a,  a'  c  L1  ( 0,“ ) j  a(t)  >  0,  a'(t)  <  0,  a"(t)  >  0,  0  <  t  <  -  (3.33) 

in  the  sense  of  measures,  and  a*  is  not  a  purely  singular  measure;  a  certain  assumption 
on  the  Laplace  transform  of  a  is  imposed  in  order  to  guarantee  that  the  third  derivatives 
of  u  are  continuous  with  values  in  L*(0,1).  The  material  function  ♦  is  also  required 
to  satisfy  4>'(0)  >  0,  and  the  technical  assumptions  regarding  the  forcing  function  f 
are  strengthened.  The  sign  of  the  memory  now  plays  a  crucial  role  in  the  local  analysis  in 
which  one  iterates  a  sequence  of  linear  integrodif ferential  equations  (compare  with  (3.18)) 
Utt  -  +  /l  a,(t-T)i).,(wj{)u)o{(x,T)dr  +  f  (3.34) 

where  u(x,t)  -  u(x,t)  for  t  <  0,  and  where  w  is  an  element  of  an  apropr lately  chosen 
function  space.  The  singular  kernel  a  satisfying  (3.33)  is  replaced  in  (3.34)  by  regular 
kernels  aj  defined  by 

a5(t)  :=■  p5  (t  )a(t+S+r  )dr  ,  0<t<-,5>0  , 

where  pj  is  a  standard  mollifler  supported  in  [-4/2, 5/2].  The  analysis 
with  singular  kernels  is  far  more  complicated  because  a"  /  L1  [0,“),  and  I  a?  I  ,  does 

IOC  o  1 

L 

not  necessarily  remain  bounded  as  540.  The  energy  estimates  are  also  considerably  more 
delicate  and  to  obtain  them  certain  technical  lemmas  concerning  Volterra  operators  with 
kernels  a  satisfying  (3.33)  are  required  (such  kernels  are  known  to  be  strongly  positive 
definite  [43]).  It  is  first  shown  that  each  linear  problem  (3.34)  has  a  unique  solution 
having  the  required  regularity  by  justifying  passage  to  the  limit  as  540.  Then  a 
contraction  mapping  argument  for  (3.34)  is  used  in  [25]  to  obtain  the  analogue  of 
Proposition  3.4  for  w  belonging  to  an  appropriate  function  space.  The  proof  in  [25]  is 
carried  out  for  B  «  [0,1]  with  Dirichlet  boundary  conditions  satisfied  at  x  -  0  and 
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x  «  1 ;  It  Is  straightforward  to  obtain  a  similar  local  result  for  B  -  R,  because  the 


local  existence  proof  in  [2S]  avoids  the  use  of  Poincard  inequalities* 

Using  their  local  result,  Hrusa  and  Renardy  then  obtain  an  analogue  of  Theorem  3.1  for 
the  history  value  problem  (2.9),  (2.10)  and  the  (singular)  kernel  a,  defined  by  (3.1), 
satisfying  (3.33)  on  bounded  Intervals.  They  impose  the  requirement  that  the  history  and 
the  solution  satisfy  Dirichlet  boundary  conditions  at  x  »  0  and  x  »  1  and  that  the 
history  and  forcing  term  be  suitably  small.  Their  result  ([25,  Theorem  5.1])  la  then  a 
simple  extension  of  the  proof  of  [13,  Theorem  1.1]  involving  the  modification  of  only  one 
estimate  in  [13];  the  modification  uses  a  refinement  of  Leona  4.2  in  [49],  because 
a"  f  [0,“)  whenever  a  is  singular.  The  fact  that  a"  ^L1[0,«>)  makes  it  difficult  to 
prove  Theorem  3.1  for  singular  kernels  using  the  analysis  in  [25].  It  is  a  challenging 
open  problem  to  prove  such  a  result  for  singular  kernels  on  all  of  space. 


4.  Development  of  Singularities  and  Related  Problems.  In  this  section  we  consider 
the  Cauchy  problem  (3.5),  (3.6)  for  regular  kernels  a,  and  we  discuss  the  development  of 
singularities  in  smooth  solutions  in  finite  time  for  smooth  but  large  data  by  using  the 
method  of  characteristics.  To  avoid  technical  complications  we  assume  that  the  forcing 
term  f  =  0  in  (3.5),  and  we  study 

ufct  -  ♦<ux)x  +  a'*i|i(ux)x,  *  «  *,  t  >  0  ,  (4.1) 

u(x,0 )  -  Ug  ( x ) ,  ut(x,0)  -  u,(x),  x  €  R  ,  (4.2) 

where  *  denotes  the  time  convolution  on  [0,tJ.  The  following  result  was  recently 
established  by  M.  Renardy  and  the  author  [42],  and  independently  by  Dafermos  [10]  for 
general  memory  functionals  using  a  somewhat  different  proof.  The  result  can  also  be 
established  by  extending  techniques  of  F.  John  [27]  to  quasillnear,  first-order  hyperbolic 
systems  with  lower  order  source  terms;  however,  the  approach  outlined  below  is  more  direct. 


Theorem  4. 1 .  Let 


€  C  («)  let  t  satisfy  (3.19)  and  let  a  be 


smooth  with  a,  a',  a"  «  LtOc^0,**‘  In  additl-on>  let  ♦  "(0)  y  0.  Then  for 

1  ?  00 

every  T ^  >  0,  there  exists  initial  data  ug,  u^  c  C  (K)  O  L  (*)  such  that  the  maximal 

Interval  of  existence  of  the  smooth  solution  u  of  the  Cauchy  problems  (4.1),  (4.2)  cannot 

exceed  T^.  More  precisely,  If  au;  |u^(x)|  and  sup  lu^Cx)!  are  sufficiently  small, 

xcl  xeS 

while  Uq ( x )  and  u^(x)  are  sufficiently  large  (with  appropriate  signs),  then  there 

* 

exists  a  number  t  <  such  that 


•up  *  t!«  <*»t)|  ♦  |u  (x,t)|)  -«  , 

l*[0,t  )  xt 


•up*  { I u  <x,t)|  +  |u  (x,t)|}  <-  I  (4.4) 

*<[0,t  ) 

For  the  special  ease  *  =  4,  Hattorl  [21]  has  shown  that  If  ♦"  ?  0  and  If  the  body 
B  is  bounded,  then  there  exist  data  u0,  u1  such  that  the  Initial-boundary  value  problem 
(consisting  of  (4.1),  (4.2)  and  compatible  Oirichlet  boundary  conditions)  does  not  have  a 
globally  defined  smooth  solution.  However,  his  method  does  not  enable  him  to  characterize 
the  data.  Ramaha  [45]  has  recently  obtained  a  blow-up  result  when  ♦  =  ♦ . 

For  first-order  model  problems  with  fading  memory,  blow-up  results  similar  to  Theorem 
4.1  have  been  obtained  by  a  number  of  authors  ([38],  [36],  [9])  by  the  method  of 
characteristics.  Existence  of  classical  solutions  for  small  data  for  such  models  is 
discussed  in  [41].  The  elegant  method  of  Dafermos  [9]  avoids  use  of  characteristics! 
Instead  a  maximum  principle  is  obtained  and  used. 

Remark  4.2.  The  reader  should  observe  that  in  Theorem  4.1  only  the  additional  hypothesis 
4"(0)  /  0  is  added  to  the  assumptions  guaranteeing  the  existence  of  a  local  smooth 
solution  of  (4.1),  (4.2)  (Proposition  3.4).  Ho  sign  information  on  the  kernel  a  is 
required.  Assumption  (3.19)  is  not  restrictive  because  it  is  shown  that  the  supremum  in 
(4.4)  is  in  fact  small. 

The  proof  of  Theorem  4. 1  generalizes  the  approach  of  Lax  [33]  using  the  method  of 
characteristics  and  generalized  Riemann  invariants.  We  transform  (4.1),  (4.2)  to  an 


mm 


(4.5) 


equivalent  first-order  system  as  follows.  Let  w  -  ux,  v  -  ut;  define 

0  J-  4(W)  -  z  ,  z  -a*  *1(1  (w)  , 

and  observe  that  a  is  the  stress-strain  functionals  (2.6).  Since  $'(•)  >  0,  equation 
(4.5)  can  be  solved  for  w,  w  »  ♦  -1(o+z)  :«  g(a,z),  and  g  is  a  smooth  function  on 
R  x  R.  As  long  as  the  solution  u  of  (4.1)/  (4.2)  remains  smooth,  (4.1),  (4.2)  is 
equivalent  to  the  system 

v.  -  o 

t  X 

«  *  c2(0,*)v  +  a’  (0)*  (g(o  ,z) )  a**#  (g(c  ,z) )  ,  ..  , 

t  x  (4.6 

z*t  -  -a*  (0)$(q(o,z) )  -  a**iMg(c,z))  , 


v(x,0)  -  u j ( x ) ,  o(x,0)  -  ♦(uj(x)),  z(x,0)  =  0  ,  (4.7) 

where  the  wave  speed  C(c,z)  i»  •  (g(a  ,z) )  )^2  is  a  smooth  function.  The  system  (4.6)  is 

hyperbolic  with  eigenvalues  C,  -C,  0.  He  define  generalized  Riemann  invariants  r,  s  by 

r  -  r(v,o,z)  v  +  4(<J,  z),  s  -  s(v,0,z)  v  -  *(o,z)i  #(0,z)  !-  f°  . 

Thus  v  »  I  the  correspondence  is  smoothly  Invertible  because  “  C-1  >  0. 

Observe  that  if  a*  a  0  in  (4.1),  z  =  0  and  g,  C  are  independent  of  z.  In  this 

situation  r  and  s  reduce  to  the  Riemann  invariants  for  the  system 

vt  -  cx,  at  -  ♦’(♦"1(o))vx 

which  can  be  transformed  to  the  quasllinear  wave  equation.  In  the  proof  r,  s,  z  are 
introduced  as  dependent  variables  and  (4.6)  is  replaced  by  an  equivalent  system  obtained  by 
differentiating  r,  s,  z  along  the  characteristics  C,  -C,  0  respectively.  One  then 
differentiates  the  quantities 


y  +  .  i.  -  i  ■  i 

X  C(9,Z) 


T  .«  V  -  — - — 

X  C(o,z) 


and  z 


along  the  C,  -C,  0  characteristics  respectively  (observe  that  if  a*  s  0,  p  -  rx, 
t  «  Xjj).  It  is  shown  (see  [42]  for  details)  that  to  leading  order  the  characteristic 
derivatives  of  /c  p,  /c  T  satisfy  a  coupled  system  of  Ricattl  equations  in  p  and  x 


with  coefficients  which  are  smooth  functions  of  r,  s,  z.  The  differential  equation  for 


zx  is  linear  in  p,  r,  zx,  and  it  is  shown  that  zx  grows  at  most  logarithmically. 
Blow-up  in  finite  time  is  established  by  showing  that  r,  s,  z  remain  in  a  neighborhood 
U  of  zero  up  to  the  blow-up  time,  if  they  are  small  initially  (i.e.  if 

sup(1v(x,0)|  +  |o(x,0) | }  is  small),  while  v'(x,0)  and  c'(x,0)  are  sufficiently  large 

* 

(with  appropriate  signs).  Moreover,  the  hypothesis  ♦"(0)  f  0  provides  upper  and  lower 
nonzero  bounds  for  the  coefficients  p2  and  T 2  in  the  Ricatti  equations  when  r,  s,  z 
are  in  t). 


Remark  4.3.  A  physical  interpretation  of  conclusions  (4.3),  (4.4),  coupled  with  examples 
of  Coleman,  Gurtin,  and  Herrera  [7] ,  is  that  the  strain  remains  bounded  but  its  first 
derivatives  become  infinite  as  t  ♦  t*.  Thus  Theorem  4.1  suggests,  but  does  not  prove, 
the  development  of  a  shock  front  in  finite  time. 


Remark  4,4.  Certain  models  for  shearing  flows  of  viscoelastic  fluids  can  be  analyzed  by 
the  technique  of  Theorem  4.1.  With  v(x,t)  denoting  the  velocity  of  the  fluid  in  simple 
shear,  Slemrod  [4B]  studies  the  problem 

vt  »  a**(vx)x»  x  £  *,  t  >  0 

(4.8) 

v(x,0)  -  v„(x)  ,  x  c  K 

in  the  special  case  a(t)  «  e-^.  Differentiation  of  the  equation  leads  to  a  Cauchy  problem 
of  the  form  (4.1),  (4.2).  Global  existence  for  smooth,  small  data  follows  from  [12, 

Theorem  4 . 1 ] >  see  also  Remark  3.3.  Development  of  singularities  for  large  data  is  an  easy 
application  of  Theorem  4.1  above.  Other  popular  models  for  viscoelastic  fluids  can  be 
discussed  by  a  similar  analysis.  Slemrod  [47]  and  Gripenberg  [20]  established  similar 
results  for  a  different  model  of  shearing  flows  for  a  viscoelastic  fluid.  If  a  -  e-t, 
(4.8)  as  well  as  the  problem  studied  in  [47],  can  be  transformed  to  the  quasilinear  wave 
equation  with  linear  frictional  damping  for  which  finite  time  blow-up  for  large  data  can  be 
established  by  the  method  of  Lax  [33]. 


Me  close  this  section  by  discussing  a  number  of  open  problems. 


Remark  4.5.  The  techniques  of  proof  of  neorem  4.1  and  that  of  [10]  depend  crucially  on 
the  hypothesis  ♦"(0)  ?  0.  The  physically  important  situation  ♦"(0)  »  0,  permitted  in 
the  finite  time  blow-up  result  for  the  quasilinear  wave  equation  (2.2]  (with  f  =  0)  in 
[37],  constitutes  an  interesting  open  problem  for  (4.1),  (4.2). 

Remark  4.6.  Singular  kernels  a  satisfying  (3.2)  and  (3.3)(ii)  -a,(0+’)  *  violate  the 
hypothesis  a"  c  L^^IO,0*)  which  is  crucial  to  the  technique  of  proof  of  Theorem  4.1  and 
that  of  the  similar  result  in  [10].  Indeed,  there  is  strong  evidence  based  on  the 
following  arguments ,  that  there  may  exist  singular  kernels  a  such  that  (4.1)  would  have 
globally  defined  smooth  solutions,  even  if  the  data  are  arbitrarily  large.  These  arguments 
suggest  that  singular  kernels  strengthen  the  dissipation  induced  by  the  memory.  Thus  far 
it  has  not  been  possible  to  resolve  this  Important  open  problem. 

First,  for  smooth  kernels  with  -a*(0+)  finite,  it  follows  from  (2.11)  and  the 
definition  of  the  constant  B  that  the  diameter  of  the  set  of  points  qQ  >  B/A  for 
which  q(t)  +  +*•  in  finite  time  shrinks  as  m(0)  ”  -a'(0+)  >0  is  increased.  However, 
the  derivation  of  (2.11)  rests  on  the  assumption  that  m(0)  «  -a*(0+)  remains  finite. 
Second,  there  are  interesting  results  of  Hrusa  and  Renardy  [26]  in  their  analysis  of  wave 
propagation  in  linear  visco-elasticity.  They  study  the  linear  history  value  problem  (2.9), 
(2.10)  with  ♦■(•)  S  c2  -  8  +  /g  m(T)dr  and  *•{•)  5  1,  u(x,t)  5  0,  t  <  0,  B  -  R,  and 
they  adjoin  step  jump  initial  data  u(x,0),  ut(x,0),  x  c  R.  They  prove  that  if  the 
memory  m  is  smooth  on  (0,-),  the  solution  has  discontinuities  propagating  along 
characteristics  of  the  linear  wave  equation  u^  “  c2uXx  an<J  *  stationary  discontinuity  of 
higher  order  at  the  initial  step- jumps.  For  singular  memory  kernels  the  propagating  waves 
are  smoothed  out.  The  degree  of  smoothing  increases  as  the  kernel  becomes  more  singular i 
the  stationary  discontinuities  remain. 

Remark  4.7.  There  is  numerical  evidence  concerning  the  development  of  singularities  in 
finite  time  for  regular  kernels  a  and  large  smooth  data.  Markowich  and  Renardy  [39]  used 
the  Lax-Wendroff  method  to  discretize  the  hyperbolic  part  in  (4.1)  and  the  trapezoidal  rule 
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to  discretize  the  integral.  They  show  that  the  method  is  second-order  convergent  and 
stable  on  any  finite  time  interval  on  which  smooth  solutions  exist.  For  spatially  periodic 
and  small  Cauchy  data,  and  for  kernels  a  which  are  finite  sums  of  decaying  exponentials, 
they  prove  second  order  convergence  on  (0,“»).  "Hiey  also  carry  out  numerical  experiments 
in  the  special  case  ^  =  4  which  exhibit  the  formation  of  a  singularity  in  finite  time  for 
particularly  chosen  ♦  ,  a,  and  suitably  large  uQ  and  u^>  Their  numerical  solution 
exhibits  but  does  not  prove  the  formation  of  shock  fronts  in  and  u^  at  the  critical 

time.  Other  numerical  schemes  merit  investigation. 

Remark  4.8.  Weak  Solutions.  Remarks  4.3  and  4.7  motivate  the  study  of  weak  solutions  for 
equations  such  as  (4.1),  (4.2)  governing  the  motion  of  materials  with  memory.  Except  for 
certain  special  situations  valid  for  steady  viscoelastic  fluid  flows  (Pipkin  [44]  and 
Greenberg  [17]),  there  is  no  rigorous  theory  for  the  existence  of  shock  waves  and 
acceleration  waves.  MacCamy  [36],  Greenberg  and  Hsiao  [19]  have  studied  several  aspects 
of  weak  solutions  but  only  for  a  single  first-order  conservation  law  with  memory  in  one 
space  dimension.  Dafermos  and  Hsiao  [11]  proved  the  existence  of  weak  solution  of  one¬ 
dimensional  first-order  quasllinear  hyperbolic  systems  with  memory  using  Glimra’s  modified 
random  choice  method  [16]  with  fractional  steps.  However,  their  method  requires 
assumptions  of  "diagonal  dominance"  which  are  not  satisfied  in  the  case  of  the  Cauchy 
problem  (4.1),  (4.2)  modelling  a  viscoelastic  solid.  They  are  satisfied  for  certain  models 
of  heat  flow  (see  [12])  and  the  specific  model  (4.8)  for  viscoelastic  fluid  flow). 

In  order  to  address  the  problem  of  weak  solutions  which  would  include  one-dimensional 
problems  for  viscoelastic  solids  of  the  form  (4.1),  (4.2),  a  program  has  been  initiated 
involving  analytical  techniques,  the  design  of  numerical  algorithms  and  numerical 
experiments.  We  consider  the  Cauchy  problem  (4.1),  (4.2)  in  the  form  of  a  first-order 
equivalent  system.  Let  w  »  ux,  v  ■  u^.  For  classical  solutions,  (4.1),  (4.2)  is 
equivalent  to  the  system 

wt  ”  vx 

vt  »  4(w)x  +  a' *i|»(w)x 


(4.9) 


satisfying  the  initial  conditions 


w(x,0)  -  wQ(x),  v(x, 0 )  -  vQ(x)  .  (4.10) 

It  is  easy  to  show  that  a  weak  solution  (in  the  sense  of  distributions)  of  (4.1), 

(4.2)  is  a  weak  solution  of  (4.9),  (4.10).  It  is  straightforward  that  the  Rankine-Hugoniot 
jump  conditions  for  elastic  shocks  (a=0  in  (4.1))  are  also  necessary  for  viscoelastic 
shocks • 

The  Riemann  problem  is  only  partially  understood  for  scalar  first-order  conservation 
laws  with  memory  (36),  but  not  at  all  for  the  viscoelastic  problem  (4.9),  (4.10). 

Therefore  it  is  difficult  to  use  the  random  choice  method  (16).  If  $  /  ♦  >  define  z  - 
a'h|>(w).  Then  (4.9)  transforms  to  the  hyperbolic  system  with  lower  order  source  terms: 

wt  -  vx 

vt  “  ♦<“>x  +  *X  (4.11) 

Zt  -  a'(0)*(w)  +  a“*(w)  , 

with  w(x,0),  v(x, 0 )  satisfying  (4.10)  and  z(x,0)  =  0.  If  ♦'(•)  >  0  (4.11)  has  the 

eigenvalues  i  (♦'(•))//2  and  0.  If  ♦'(•)  -  a(0)4'*(“)  >  0,  (4.11)  has  a  uniquely 

determined  steady  state  solution.  Observe  that  initially  *x  =  0i  one  can  solve  the  first 

two  equations  in  (4.11)  by  various  techniques  for  conservation  laws  on  the  first  time  step, 

update  z  using  the  last  equation  and  proceed  forward  in  time.  Jointly  with  B.  Plohr  we 

have  initiated  a  study  of  various  numerical  algorithms  for  (4.11)  in  the  special  case 
n 

a(t)  “  1.  \  exp(-X  t),  a^  >  0,  1^  >  0,  including  the  Glimm  scheme  with  fractional 

i»1  K 

steps.  One  objective  is  to  establish  existence  of  weak  solutions  for  small  BV  data. 
Another  is  to  obtain  implementable  numerical  algorithms  which  can  be  tested  on  concrete 
problems . 

Boldrlni  (3,  4]  used  techniques  of  compensated  compactness  to  study  elastic  and 
viscoelastic  problems  including  the  system  (4.9),  (4.10).  these  techniques  were  developed 
by  Tartar  (50,51,52),  Murat  (40)  and  DiPerns  (14) t  in  [14]  DiPerna  succeeded  to  extend 
these  techniques  and  apply  them  to  establish  the  existence  of  weak  solutions  of  the  purely 
elastic  one-dimensional  problem  (i.e.  (4.9),  (4.10)  with  a  =  0)  on  Rx  [ 0 , T)  for  any 
T  >  0,  without  restricting  the  size  of  the  data.  Boldrini  [4J  assumes  that  the  memory  in 
(4.9)  is  small  in  the  sense  that 

a  a(5,t),  ♦(•)  ♦(•)  +  ug(*  )  ,  (4.12) 

where  4  >  0,  u  >  0  are  small  parameters,  g  is  a  smooth  function  satisfying  the  growth 
condition  |g(w)|  <  x|w|  ,  k  >  0,  and  a'(5,t)  =  0(5),  a"(6,t)  »  0(5)  uniformly  in  t. 

In  place  of  (4.9)  he  considers  the  regularized  system 


vt  »  4(w)x  +  a*(5/*)  *  (4(w)  +  ug(w))x  +  evxx 


(4. 13) 


with  initial  data  (4.10)  (the  Newtonian  viscosity  can  be  more  general  than  Evxx)»  where 

£  >  0  is  a  small  parameter.  Let  w  ,  ,  v  ,  be  a  solution  of  (4.9),  (4.10)  on 

e ,5 ,u  £ ,6  ,u 

R  »  [0,TJ  for  any  T  >  0.  Boldrini  gives  sufficient  conditions  which  insure  that  there  is 

a  subsequence  such  that  w  ,  ♦  w,  v  ,  -*■  v  on  R  x  [0,T]  as  e  ,  5  ,  u  *  0+,  where 

£ ,o  ,U  £ ,U 

w  *  Moreover,  w,  v  is  a  weak  solution  of  the  purely  elastic  problem  on 

R  x  [o,T] .  The  most  serious  of  his  assumptions  is  the  crucial  hypothesis  requiring  the 
solutions  w.  t  v,  i  ,,  of  (4.13)  to  lie  in  l“  uniformly  in  the  parameter  £,  6,  u. 
Since  the  memory  i3  a  nonlocal  operator,  this  assumption  is  difficult  to  verify. 

Jointly  with  W.  Rogers  and  T.  Tzavaras,  we  are  using  compensated  compactness 
techniques  to  establish  the  existence  of  weak  solutions  of  (4.9),  (4.10).  The  special 
case  p  =  p,  but  with  the  memory  not  small  (i.e.  a  independent  of  6  )  is  tractable  by 
these  methods  and  the  case  <|i  f  p  appears  doable.  However,  obtaining  an  invariant  region 
in  order  to  show  that  solutions  of  the  relevant  regularized  system  lie  in  l“  is  extremely 
difficult.  It  is  of  interest  to  note  that  the  existence  of  weak  solutions  of  the  Cauchy 
problem  for  the  model  first-order  scalar  equation  with  memory 
ut  +  *<u)x  +  a'  *i(i  (u)x  -  0,  x  £  R,  t  >  0 


(4.  14) 


U ( X, 0 )  “  Uq ( X )  ,  X  €  R  , 

where  a,  P,  P  have  the  same  meaning  as  in  (4.9),  can  be  solved  completely  by  using  the 
method  of  compensated  compactness.  The  maximum  principle  proved  by  Dafermos  in  [9]  for 
classical  solution  of  (4.14)  makes  it  possible  to  prove  the  needed  l“  estimates  for 
solutions  of  the  regularized  problem  (i.e.  (4.14)  with  Euxx  on  the  right  side  in  place  of 
zero).  This  problem  was  recently  solved  by  Dafermos  (oral  communication).  Unfortunately, 
it  does  not  appear  that  this  approach  can  be  extended  to  coupled  two  by  two  systems  with 
memory. 
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